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Abstract 

A series of tests has been performed on a four-port wave rotor suitable for use as a topping stage on a gas 
turbine engine, to measure the overall pressure ratio obtainable as a function of temperature ratio, inlet mass flow, 
loop flow ratio, and rotor speed. The wave rotor employed an open high pressure loop that is the high pressure inlet 
flow was not the air exhausted from the high pressure outlet, but was obtained from a separate heated source, 
although the mass flow rates of the two flows were balanced. This permitted the choice of a range of loop-flow 
ratios (i.e., ratio of high pressure flow to low pressure flow), as well as the possibility of examining the effect of 
mass flow imbalance. Imbalance could occur as a result of leakage or deliberate bleeding for cooling air. 
Measurements of the pressure drop in the high pressure loop were also obtained. A pressure ratio of 1.17 was 
obtained at a temperature ratio of 2.0, with an inlet mass flow of 0.6 Ib/s. Earlier tests had given a pressure ratio of 
less than 1.12. The improvement was due to improved sealing between the high pressure and low pressure loops, 
and a modification to the movable end- wall which is provided to allow for rotor expansion. 

Nomenclature 

G leakage parameter = 25/h 

H passage height 

A:ioss loss coefficient for heater loop 

m mass flow rate 

N rotor speed (rpm) 

P stagnation pressure 

T stagnation temperature 

TR target temperature ratio for TJTx 
U mixed-out value of velocity 

y ratio of specific heats 

AT temperature change due to heat transfer 

5 gap between rotor and end- wall 

T| efficiency 

X loop flow ratio (= 

p mixed-out value of density 
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Subscripts 
c compression 

X expansion 

1 Port 1 

2 Port 2 

3 Port 3 

4 Port 4 


I. Introduction 

At NASA, wave rotors have been investigated for use as a topping stage for gas-turbine engines in order to 
increase performance while maintaining temperatures in rotating parts at levels used in current technology. The 
primary application considered is aircraft propulsion or auxiliary power. Calculations have shown that increases in 
specific power of ~20 percent, and reductions in specific fuel consumption of ~15 percent are possible (refs. 1 and 
2). The original intended use of the wave rotor was as a topping stage for ground-based power units (ref 3) and it 
has been demonstrated that large increases in efficiency and power are possible for that application (Zauner et al.) 
(ref 4). In the configuration used by Zauner et al., the power generated by the topping stage was extracted by an 
additional power turbine, with its own shaft. Thus this power unit with a wave rotor topping stage is more complex 
than a power unit without the topping stage. 

For aircraft engines, a simpler, preferable, arrangement would be to have the extra power be produced by an 
increased pressure at the turbine entrance. This arrangement is schematically illustrated in figure 1 . In this scheme, 
the wave rotor acts on a four-port cycle. Air from the engine compressor enters the wave rotor at port 1, and its 
pressure is further increased inside the wave rotor by unsteady compression waves. The compressed air leaves the 
wave rotor at port 2, and passes to the burner. Combustion therefore takes place at higher pressure than in the engine 
without a wave rotor. From the burner, the gas re-enters the wave rotor at port 3, expands in expansion waves, and 
leaves the wave rotor at port 4, passing to the engine turbine. This gas is now at higher pressure than the gas entering 
the turbine in an engine without a wave rotor, even though the temperature is the same, due to the cooling in the 
expansion. Consequently, more work can be extracted from the gas in the engine with the wave rotor. 

An investigation of an application of a wave rotor topping stage to an existing aircraft engine, using a four-port 
wave rotor cycle has been reported by Snyder and Fish (ref 5). The performance of the engine with the wave rotor 
topping stage depends on the performance of the wave rotor itself, which in this context means the pressure ratio 
developed across the wave rotor as a function of the temperature ratio across it, i.e., P^IPi versus T^/Ti. 

There appear to be only two investigations of four-port wave rotor performance, by Klapproth et al., (ref 6) and 
by Ruffles (ref 7) (see figure 2, which also includes a one-dimensional calculation of wave rotor performance 
(ref 8). There is significant discrepancy between these two reports, and neither is very complete. Thus there is a 
need for detailed measurements of four-port wave rotor performance. Such measurements are also needed to 
validate numerical simulations. This paper gives the results of an experiment to measure four-port wave rotor 
performance. 


Burner 



Figure 1. — Schematic arrangement of a wave rotor in a gas turbine engine. 
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Figure 2. — Experimental results of Klapproth et al., and Ruffles, together with 
predictions from a 1-D code calculation by Paxson (ref. 8). 
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II. The Wave Rotor Cycle 

As stated above, the 
compression and expansion in 
a wave rotor are performed by 
waves. The rotor itself is a 
cylindrical drum, with many 
axial passages on its 
circumference, and the waves 
propagate within these 
passages. A convenient way to 
represent these processes is a 
so-called x-t diagram, which 
is a diagram showing the 
position of the waves as a 
function of time, i.e., as the 
cylinder rotates. An example 
is given in figure 3, showing 
the temperature (fig. 3(a)) 
and pressure (fig. 3(b)), as 
calculated with the one- 
dimensional wave rotor code (see ref 8, and the appendix for details), for one of the runs in the test series. The 
temperature and pressure are plotted along a passage of the wave rotor as it rotates around, passing the various ports, 
and returning to its original position after one complete revolution. The relative width of a passage is indicated at the 
top of the diagram. Also shown, as grey circles, are the probes in port 4. Starting at the bottom, the passage is fairly 
uniform in temperature, except for a small region of cold gas at the right hand end. The pressure not quite uniform as 
a shock wave is propagating from left to right, and is approximately two-thirds of the way along at time zero. At 
time zero, the passage passes the leading edge of port 4 at the right hand end. Port 4 is at lower pressure than the gas 
in the passage, so that an expansion wave propagates into the passage, traveling to the left, and reducing the 
temperature. When this wave reaches the left hand side is when the passage has rotated until its left hand end is at 
the leading edge of the inlet port. The pressure in the passage has been reduced by the expansion wave, and is lower 




Figure 3.- 


-The space-time development of waves inside the wave rotor, 
(a) Temperature, (b) Pressure. 
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than that in port 1 , so cold inlet air flows in. As the passage continues to rotate, the right hand end passes the trailing 
edge of port 4, and so is closed off. This generates a compression wave (hammer shock) to bring the flow to rest, 
which propagates to the left, reaching the left hand end as the passage has rotated to the trailing edge of port 1 , 
trapping high pressure air in the passage, since it is closed at both ends. Next the passage rotates until the left hand 
end is open to port 3. Port 3 is at high pressure and temperature, having come from the heater, so flow enters the 
passage, sending a shock wave to the right, compressing the air in the passage even more. When this shock wave 
reaches the right hand edge, the passage has almost rotated to the leading edge of port 2. The shock wave is reflected 
from the wall, propagating back to the right, but is followed by an expansion wave, as the high pressure gas escapes 
into port 2. The air in port 2 flows to the heater. Most of the air that entered from port 1 flows into port 2, but a 
portion has not flowed out when the right hand end of the passage is closed by reaching the trailing edge of port 2. 
Earlier, when the left hand end of the passage reached the trailing edge of port 3, this caused the flow to stop, 
creating an expansion wave, propagating to the right. When this wave reaches the right hand end of the passage is 
the time at which the passage has reached the trailing edge of port 2. The passage is now closed off again, and 
reaches almost uniform conditions (except for the cold air trapped at the extreme right hand end) before the whole 
process starts again. This is the so-called through-flow cycle, which has the advantage that the rotor is exposed 
alternatively to hot gas and cold air, and thus is self cooled. Other cycles are possible, but this cycle was chosen as 
appropriate for the application under consideration. 

III. Description of the Wave Rotor 

The heart of a wave rotor is the rotor itself. The rotor used initially is shown in figure 4. It was 8 in. outer 
diameter, 10 in. long, with webs 0.030 in. thick. The dimensions of the rotor were determined using an optimization 
procedure (ref. 1), and the one-dimensional wave rotor computer code developed by Paxson (ref 8). The rotor was 
constructed from titanium. The 
passages were first drilled, then 
electric-discharge machined to 
final shape. The web design 
proved disastrous, as it failed 
after about 2 hr of running. A 
second rotor was built with more 
robust webs, which were also 
chemically etched to remove the 
recast surface layer, and the root 
areas shot-peened to put them in 
compression (ref 9). The webs 
can be seen in figure 5, which 
shows the rotor with the movable 
end- wall on top of it. A movable 
end- wall is provided between the 
rotor and the fixed end-wall at 
the outlet end. The purpose of 
this feature was to provide for 
thermal expansion of the rotor. 

The end-wall had a fairly tight 
clearance between it and the 
rotor, and was located in a fixed 
position on the rotor shaft Figure 4. — The original NASA wave rotor, 

relative to the rotor at the outlet end. The rotor was 

fixed relative to the casing at the inlet end, so all thermal growth occurred at the outlet end. A fairly large gap was 
designed between the movable and fixed end-walls at the outlet end to allow for this expansion. This gap was sealed 
axially by a moustache seal on its rim, and circumferentially, i.e., to prevent flow from high pressure to low pressure 
inside the gap, by spring vanes on the fixed end-wall. In addition, graphite seals were installed to prevent radial flow 
from the rotor to the cavity around the rotor. These can not be seen in the cross-sectional view of the rig given in 
figure 6, but were placed in the open squares seen at each end of the casing outside the rotor. 
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Figure 5. — View of movable end-wall sitting on the rotor. 



Figure 6. — Cross-sectional view of the assembled wave rotor rig. 


The rotor is housed inside a casing, between the end- walls, to which the port ducts are attached. A photograph 
of the assembled rig is given in figure 7. The ducts bring the air in or out depending on the particular port, and are 
set at the angle such that the circumferential component of the gas velocity is equal to the velocity of the rotor at the 
radius corresponding to the centerline of the passages. A photograph of the port 1, inlet, duct is given in figure 8. 
The ducts are transitions between the port shape, and the circular shape of the piping. The exit ducts, which are 
simultaneously diffusers, are necessarily much longer. 

The rotor is driven by a variable speed 10 hp electric motor, governed by a constant speed controller, with a 
maximum speed of 12,500 rpm. 
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Figure 7. — Photograph of the assembled wave rotor rig. 



Figure 8. — Photograph of the wave rotor inlet duct. 


IV. The Experiment 

In figure 1, the high pressure exhaust (leg 2), is fed to the burner, and from there back into the wave rotor at 
port 3. Although the experiment was initially set up this way, it did not work, as the heater pressure losses were too 
large. Instead, the experiment was run as an open loop. That is to say that the exhaust at port 2, after passing through 
the diagnostic spool, and the orifice meter, was simply dumped to vacuum. A supply of high pressure air, after being 
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heated, was sent through the venturi meter, and to port 3. This is shown in figure 9. Both legs were equipped with a 
valve to control the flow rate. In running the experiment, the operator adjusted the valves until the mass flow rates 
were equal, and matched the value requested for any given run. This gave the ability to choose different values of X. 
It was found that this ratio is related to the pressure drop between port 2 and 3. Thus, being able to vary X proved to 
be an advantage as it permitted exploration of the effect of differing values of loop loss coefficient on wave rotor 
performance. The loop loss coefficient is defined in the normal manner as: 

hoss={P2-P3)X-5p2U^ ( 1 ) 

The experiment differed from figure 1 in another respect. Rather than a burner, a Caloritech T90-2444, 350 kW 
electric heater was available, and was used. Use of a heater has an advantage and a disadvantage. The advantage is 
that the flow corresponds to currently available computer codes for calculating wave rotor performance (ref 8), 
which are limited to a single value of the ratio of specific heats. This is simultaneously its disadvantage, in that a 
burner would produce a gas with a ratio of specific heats more appropriate to that produced in an engine, and 
therefore more true to life. Another disadvantage became apparent after a few runs. The heater was designed for a 
different experiment, involving similar mass flows, but at much higher pressures. When operated at the pressure of 
the present experiment, the stagnation pressure drop across the heater was very large, which led to the decision to 
use an open loop. 

The experiment was run entirely at an input pressure T’l of 7.5 psi, and an inlet temperature of 400 R. The latter 
value was used to limit the temperatures necessary in the heater, while still providing temperature ratios up to 
T 4 /T 1 = 2.0. Different schemes for achieving this were considered (ref 10), with the present arrangement being 
selected based on the availability of a chilled air supply at NASA Glenn Research Center. The pressure was chosen 
to reduce the forces on the rotor webs which had led to the failure of the original rotor, which was run at an inlet 
pressure of one atmosphere. 


V. Instrumentation 

Figure 9 is a schematic view of the experiment. A diagnostic spool is installed at the outer end of each duct, 
with three static pressure taps, plus a rake with five total pressure probes (each at a different radial location), and 
two thermocouples. From these pressure and temperature measurements, the radial velocity distribution can be 
found, and hence an average total pressure, at each duct. Similar ducts were provided at the entrance to, and exit 
from, the heater. Tube-type combination pressure probes (ref 1 1) are used in the inlet port to measure flow direction 
and, in conjunction with wall static taps and the known inlet temperature, velocity. In the other ports, rakes with 
three pitot probes, and two thermocouples, together with wall static taps, are used to give radial velocity and 
temperature distributions. Four such rakes are used in the exit port 4, plus another rake with five pitot probes to give 
a more complete radial velocity distribution. The resulting circumferential and radial velocity distributions are used 
in a mixing calculation (ref 12) to give the exit average total pressure and temperature. All pressure measurements 
are recorded through an electronically scanned pressure (ESP) measurement system. The pressure measuring system 
automatically self-calibrates every hour to maintain a 0.1 percent accuracy. In the exit leg, a temperature spool, with 
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nine thermocouples, each sampling an equal section of the spool area, with the resulting signals electronically 
averaged, is used to give the exit temperature T^. Thermocouple readings are corrected to total temperature readings 
using the recovery correction factor of Glawe, Krause, and Dudzinski (ref. 13). The rotor speed is read by a once- 
per-revolution proximity counter. All the above signals are fed to a computer, which takes 30 readings for each data 
point, and then stores the average value. 

Mass flow rates are measured by standard orifice meters on the inlet leg 1, and the outlet legs 2 and 4. The mass 
flow rate from the heater is read by a venturi meter. Since venturi meters are more accurate than orifice meters, each 
orifice meter was calibrated against the venturi meter, by arranging for flow to pass only through the venturi meter 
and the orifice meter in calibration. The flows from ports 2 and 4 are very non-uniform in temperature. Each orifice 
meter uses one thermocouple to monitor the temperature of the air passing through it. It is possible then that the 
mass flow rate value could vary depending on whether the thermocouple sampled the hot portion of the flow, or the 
cold portion. After some early runs showed considerable variation in results, a VORTAB© was installed upstream 
of the leg 2 orifice meter, to provide significant mixing of the flow, thereby ensuring a single value of temperature. 
This appeared to solve the problem. 


VI. Data Corrections 

Despite the best efforts of the operators, there was inevitably some imbalance in the high pressure loop mass 
flows, and also for those runs in which a constant temperature ratio was required, some variation in temperature 
ratio. Consequently, for any given series of readings, the data was entered into a statistical analysis which included 
(w3 /m2 ~l) ? and {T/[ITi-TR) as error parameters in addition to the parameters affecting the change in performance 
being investigated. If the variation due to these errors was significant, the program would provide a multiplying 
constant, which was then used to correct the data to give a value corresponding to zero error. The variation in rotor 
speed was so small that there was no need to correct for any variation. 

VII. Results 

The test series consisted of three basic types of tests. These were (a) measurement of the ratio Pa/P\ at a fixed 
inlet mass flow as the temperature ratio was raised, (b) measurement of the ratio P^/Pi at a fixed value of r4/ri, and 
fixed values of N and X, as m\ was changed, for different values of X, and (c) measurement of P^/Pi at fixed T^/Ti, 
m\ and X, as N was changed. 

A. Effect of Variation of Temperature Ratio 

Runs were made in which the pressure ratio P^/Pi was measured as TJTi was raised from about 1.6 to 2.0, at 
various values of m\^ m 2 = m^ =1.08 Ib/s , and a fixed A = 10,600 rpm. This is not representative of a real wave rotor 

with port 2 connected directly to port 3. Such a wave rotor would have a pressure drop in the upper loop determined 
by the pressure loss factor, i.e., A:ioss, plus a pressure drop due to heat addition (ref 14). The latter is given 
approximately by: 

{P^-P2)hoth■^P2Ul={T^IT2-\) (2) 

and is typically small compared with the pressure drop due to the pressure loss factor. Running the experiment at 
constant A:ioss would therefore correspond roughly to a real wave rotor. This was in fact the procedure for some very 
early runs made with the present wave rotor. However, it was difficult to do, since this involved setting a different 
value of X for each point as the temperature was raised. These runs had shown that the performance did not vary 
much with variation in X, and so it was decided to fix the value of m 2 = ^3 in the present tests, as this made the 

tests much easier to run. The value chosen, 1.08 Ib/s, corresponds to the design value of X = 1.6 for an inlet mass 
flow of 0.67 Ib/s. The results of the tests are given in figure 10. The tests were run for a range of inlet mass flows, 
from the maximum possible of 0.72 Ib/s when the valve V4 is wide open, to the lowest value of 0.55 Ib/s. The lines 
are least square fits to the data points, which are all of the form: 

P3 !P\=a + b{T4 /Ti )+c{T 4 /n f (3) 

The values of a, b, and c for each inlet mass flow rate are given in table 1. Also given in figure 10 is a curve of 
pressure ratio versus temperature ratio for an inlet mass flow of 0.69 Ib/s calculated with the one-dimensional wave 
rotor code, showing remarkably good agreement. 
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74/71 

Figure 10. — PJP\ versus TJTi at various inlet mass flow rates. 


TABLE 1.— VALUES OF CONSTANTS IN EQUATION (3) 


m\ 

0.72 

0.69 

0.67 

0.63 

0.59 

0.55 

a 

-0.358 

-1.662 

-1.270 

-0.663 

0 

0 

b 

1.192 

2.645 

2.326 

1.768 

1.041 

1.03 

c 

-0.233 

-0.618 

-0.554 

-0.426 

-0.226 

-0.221 


B. Effect of Variation of Inlet Mass Flow Rate 

The inlet mass flow rate was varied for a fixed rotor speed, V= 10,600 rpm, at TR = 2.0, for each of three 
values of ^2 = ^ 3 , namely 0.9, 1.08, and 1.2 Ib/s. The results are given in figure 11. Fits to the points are also 
given. The performance increases rapidly at first as the mass flow is reduced from the maximum possible, then tails 
off after the flow rate reaches 0.67 Ib/s. The lowest value of m 2 used provides a lower performance. Lines of the 
pressure drop in the high pressure loop, i.e., {P 2 -P^)IP 2 are also plotted in figure 11. The pressure drop increases 
rapidly as m 2 decreases. It might be thought that a low drop is desired, but the situation is more complicated — to 
obtain a high value of m 2 , a low pressure drop is required, so that the high pressure loop must have a low value of 
A:ioss- In figure 12, X (closed symbols) and {P 2 ~P 2 )IP 2 (open circles) are plotted against Not surprisingly, 
{P 2 ~P^)IP 2 increases monotonically with A:ioss, whereas X decreases. Also indicated is the minimum combustor 
pressure drop given by Lefebvre (ref. 14). For the points taken at m 2 =1.2 Ib/s, is low, and the low pressure 
drop required to achieve this is lower than the Lefebvre minimum loss. Thus this is not a realistic situation for a 
practical wave rotor. The points taken with m 2 = 0.9 Ib/s, on the other hand, have a large A:ioss, and resulting large 
pressure drop, so can be a realistic configuration, but suffer from lower performance. For a practical wave rotor, the 
burner loop should be designed to have as low a pressure loss as possible. In fact it is also necessary for the duct 
from port 4 to the engine to have low loss. Studies of a potential duct design are given in reference 15. 
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C. Effect of Variation of Rotor Speed 

Tests were run in which performance was measured as the rotor speed was varied at TO = 2.0, for each of the 
inlet mass flow rates used in subsection (B), with m 2 = m^= 1.08 Ib/s. The results, given in figure 13, show bell- 
shaped curves, with the speed of maximum performance increasing as the inlet mass flow drops, from 10,200 rpm 
for the maximum inlet flow rate, to 10,900 for the lowest flow rate. Near the peak, however, all the curves are fairly 
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Figure 13. — P 4 /P\ versus N at T 4 /T 1 = 2.0, for various inlet mass flow rates. 


flat, so that using a value of 10,600 rpm for all inlet mass flows, as was done in subsection (B), does not result in 
much loss of performance. On the other hand, a speed below 10,000 rpm, or above 12,000 rpm, does cause a loss of 
performance, particularly for the lower inlet mass flow rates. 


D. Compression Ratio and Efficiency 

It is of interest to see what compression ratio, i.e., P 2 /P\, is achieved in the wave rotor. This can be determined 
from the results of subsection VII (B), and is plotted in figure 14, showing the compression ratio for the different 
inlet mass flow rates, for each of the three values of m 2 = . The compression ratio is seen to be between 2.75 

and 3.2 depending on the inlet mass flow rate and the value of m 2 = m^ . In figure 15, the corresponding 
temperature ratio, i.e., T 2 IT 1 is plotted. If the adiabatic compression efficiency is calculated as: 

11, = ri ((P 2 jPx - i)/(r2 - ri ) (4) 


a very low value, around 30 percent is obtained. Similarly, if the expansion efficiency (from state 3 to state 4) is 
calculated, a value greater than unity is found. This is obviously incorrect. The explanation lies in the fact that the 
process is not adiabatic. The hot gas entering at state 3 contacts the rotor, which is approximately at temperature T 2 , 
and so is cooled while it is expanding, and the cold air entering at state 1 is heated by the rotor while it is being 
compressed. Thus the temperature rise (r 2 -ri) is partly due to work of compression, and partly to heat transfer. If it 
is assumed that ATc is the temperature rise due to heat transfer to the cold entering air, and ATx the temperature drop 
due to heat transfer from the hot entering gas, then the compression efficiency is really 


11c = 7l 




and the expansion efficiency is; 


= [73 -74 -Arj/73 


(5) 


( 6 ) 


once the rotor has reached a steady state temperature, there is no net heat transfer to it, and so, since the mass flow 
rate of hot air is larger than that of cold air. 


XAT^=ATc 


(7) 
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Figure 14. — Compression ratio versus inlet mass flow rate at TR = 2.0, A^= 10 600 rpm. 
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Figure 15. — T 2 IT 1 versus inlet mass flow rate at TR = 2.0, N= 10 600 rpm. 

without knowing the values of AFc, or AF^, the efficiencies can not be calculated. However a rough idea of their 
values can be found by assuming they are equal. With this assumption, equations (5), (6), and (7) can be solved to 
give Tjc = Tjx- The results are plotted in figure 16 for the tests of subsection VII (B), and the efficiency is seen to be 
around 0.66 to 0.74. This is in line with values for wave rotors given by Thayer and Zumdieck (ref 16) (0.75), 
Ruffles (ref 7) (0.77), and Kollbrunner (ref 17) (0.65 to 0.68). Kollbrunner’s result was for a direct measurement of 
compression efficiency. Values of ATx range from 200 to 300 R. 
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Figure 16. — Compression/expansion efficiency versus inlet mass flow rate at 
TR = 2.^,N= 10 600 rpm. 


However, this still does not represent the compression efficiency due to waves alone. As was pointed out in 
section II, and shown in figure 3, not all the air entering at port 1 exits at port 2; a portion of it is trapped on the rotor 
and exits at port 4. Similarly the flow into port 2 has a mixture of air from port 1 and air from port 3. The values of 
P 2 , and P 4 are thus mixed-out values of the two contributions to the port flows. The mixing inevitably also involves 
loss, which is not separated out in the above calculation of compression and expansion efficiency. 


E. //i2 ///13 Imbalance T ests 

In an actual wave rotor, the high pressure flow out of the wave rotor at port 2 may not equal the high pressure 
flow into the wave rotor at port 3, either due to leakage, or due to deliberate bleed at port 2 for cooling purposes. It is 
of interest to see how this might affect performance. A test was run in which rh^ was systematically reduced below 

m 2 and the performance measured. The test was run at the maximum inlet mass flow rate, TR = 2.0, and for each of 
m 2 = 0.9, 1.08, and 1.2 Ib/s. The results are given in figure 17. As expected, the performance is reduced when 
m 3 < m 2 The solid lines are least square fits to the data, while the dashed lines are approximate fits given by; 


Pa!P\ = A/^l(w3/w2=l) + 0.639X(m3/m2 -l) 


( 8 ) 


which seem to fit almost as well as the least squares fits. 

VIII. Discussion 

Reference was made above to earlier tests performed in July 1999. In figure 18, values of PaIP\ for an input 
mass flow rate of 0.6 Ib/s from those tests have been plotted versus T 4 IT 1 . Also shown are the results of Klapproth 
et al., (ref 6 ) Ruffles (ref 7), and the present tests, at an input mass flow rate of 0.59 Ib/s (April 2006). The 1999 
tests showed lower performance than that of Klapproth et al., and the code calculation, although better than that of 
Ruffles. The current tests are clearly superior in performance. There are several contributors to this improved 
performance. The first is reduced clearance gaps. Based on the conclusion by Wilson (ref 18) that a leakage 
parameter G < 0.010 will provide essentially leakage free performance, the gap between rotor and end-wall at the 
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m^/i7i2 

Figure 17. — The effect on performance of imbalance of mi, I m 2 , TR = 2.0, ml = 0.12. 



Temperature ratio 74/T'i 

Figure 18. — P 4 /P 1 versus T 4 /T 1 from the NASA tests (July 1999, and April 2006), 
together with results from Klapproth, Ruffles, and code calculations. 

inlet end had been set to 0.010 in. (G = 0.008) and 0.020 in. (G = 0.015) between the rotor and the movable end- wall 
at the outlet end for the early tests. These were considered to be sufficiently small. However, the gaps were reduced 
to 0.004 in. at the inlet end, and 0.002 in. at the outlet end for the later tests. Secondly, there was evidence from tests 
that the movable end- wall was tilting due to the asymmetrical pressure loading on its rear face. To reduce this tilting, 
the movable end-wall was provided with a double bearing set for the 2006 tests. Thirdly, it was realized that there 
was a considerable path for circumferential leakage in the space around the moustache seal. This was sealed with 
RTV glue in the last tests. Finally, seals were incorporated to prevent axial leakage into the cavity around the rotor. 
In the 3-port wave rotor tests reported by Wilson (ref 18), similar axial sealing was provided by brush seals, and 
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proved very effective. It was not possible to obtain brush seals for the four port wave rotor, so graphite seals were 
used instead. These were not available for the 1999 tests, but were in place in the 2006 test series. 

Although the movable end-wall seemed a good solution to the problem of differential expansion between rotor 
and casing, in practice it proved troublesome, as explained above. A better solution might be to not use a movable 
end- wall, but to design a cold clearance equal to the anticipated difference in expansion. If the expansion proves to 
be greater than anticipated, abradable material on the end- wall can provide the extra space needed. With this 
technique very tight clearances may prove possible. 

The results of subsection VII(B) indicate that wave rotor performance is linked to the losses in the high pressure 
loop. A cycle designed on wave considerations alone implies a loss in this loop, and the loop must equal this loss or 
the desired conditions will not be achieved. 


IX. Conclusions 

Wave rotor performance has been shown to be determined by many factors, namely temperature ratio across the 
wave rotor, loop flow ratio, inlet mass flow rate, rotor speed, and also leakage paths. Loop flow ratio also affects the 
pressure loss between ports two and three, and the operating point may be constrained by the loss in this leg. It is 
important to reduce the leakage as much as possible. The observed pressure ratio of 1.17 at T/^ITi = 2.0 is very good, 
although it might still be possible to get further improvement if leakage could be reduced more. Higher pressure 
ratios can obviously be achieved at higher values of 

Appendix 

The CFD code used in this paper to simulate the wave rotor has been documented extensively in the literature 
(see for example reference 8), and only a brief description will be provided here. It is a time accurate, quasi-one- 
dimensional code which numerically integrates the equations of mass, momentum, and energy for a calorically 
perfect gas, in a single passage of the wave rotor, as it rotates past the various ports governing the cycle. A second 
order, Lax-Wendroff integration scheme is utilized which incorporates Roe’s approximate Reimann solver and a 
flux limiter in order to ensure total variation diminishing behavior, while preserving sharp resolution of shock waves 
and contact discontinuities. Source terms have been added to account for the effects of leakage between the passage 
and the casing, wall friction, and wall heat transfer. The sub-models which comprise the source terms have been 
developed using results from previous wave rotor experiments. They contain correlations which link their strength to 
physical aspects of the experimental rig (e.g., rotor-to-endwall clearances, passage hydraulic diameter, Prandtl 
number, etc.), and local fluid dynamics. Source terms are applied globally for wall friction and heat transfer, and 
discretely (i.e., the first and last interior numerical cell) for leakage. 

A robust, characteristics-based boundary condition algorithm has been developed which permits inflow or 
outflow at any port, depending on the status of the passage flow and the specified port pressures. This allows 
simulation of both on and off-design operation. Inflow losses due to incidence are accounted for using sub-models. 
The code operates in the passage relative frame, while boundary conditions are specified in the absolute frame. The 
boundary conditions routines account for this reference frame change as well. 

In order to obtain meaningful limit cycle results with which to compare experimental results, two components 
have been added to the code. The first is a simple cavity representing the casing volume not occupied by the rotor. 
Leakage flow travels to and from this volume (via the code source terms) as the passage rotates. The amount and 
direction of leakage depends on the state of the cavity, and the state of the passage cells where leakage is assumed to 
occur. The gas state in the cavity is updated each revolution, or cycle, of the rotor. When the net flow to the cavity 
over the course of a cycle is zero, a limit cycle has been reached. 

The second additional component is a burner model, which links ports 2 and 3. It is again modeled as a simple 
cavity; however, there is a heat addition term and a total pressure loss term added. The cavity state is governed by 
the computed flows from port 2, and to port 3, as well as by the specified heat addition rate. The state is updated 
each cycle via first order time integration. When the flow into the burner matches the flow out, and the burner state 
is unchanged from one cycle to the next, a limit cycle has been reached. 

For the results presented in this paper, endwall-to-rotor clearances are assumed to be 0.004 in. at the inlet end 
and 0.002 in. at the exhaust end. Boundary conditions imposed are as follows. The measured stagnation pressure and 
temperature are specified for port 1. Average static pressure is imposed for port 4. It is adjusted until the computed 
inlet flow matches the experimentally measured value plus approximately 8 percent to account for blockage due to 
the passage ‘webs’. These are not modeled in the code. Total conditions for port 3 and static pressure for port 2 are 
supplied by the burner component model. The heat addition term and the total pressure loss terms are adjusted until 
the computed and measured values of TR and X match. 
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